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Applications  of  Sphere  Geometry 
in  Canal  Surface  Design 


Christoph  Maurer 


Abstract.  Classical  models  of  sphere  geometry  facilitate  an  efficient 
description  of  canal  surfaces.  Using  the  cyclographic  model  of  Laguerre 
geometry,  an  elementary  characterization  of  continuity,  bisectors,  control 
point-,  control  sphere-  and  implicit  representations  of  canal  surfaces  is 
presented.  In  addition,  canal  surfaces  generated  with  the  aid  of  Minkowski 
Pythagorean  hodograph  curves  are  investigated. 


§1.  Introduction 

A canal  surface  C in  IR^  is  defined  as  the  envelope  surface  of  a moving  sphere 
S{t)  with  center  m{t)  :=  {mi{t),m2{t),m^{t))  and  radius  function  r{t).  The 
moving  sphere  can  be  described  with  the  implicit  equation 

F(*,t)  = l|®-m(t)|p-r2{t)=0.  (1) 

The  envelope  condition 

dF 

- = [x-  m(t)]  • rh  + r{t)r{t)  = 0 (2) 

describes  a moving  plane,  which  intersects  S{t)  in  the  characteristic  circles  of 
the  canal  surface.  It  is  a natural  approach  to  use  models  of  classical  sphere 
geometry  to  study  canal  surfaces.  Well-known  models  have  been  investigated 
in  Mdbius  geometry,  Laguerre  geometry  and  Lie  geometry.  For  an  overview 
on  sphere  geometry,  the  reader  may  consult  [2].  Papers  which  handle  sphere 
geometry  in  the  CAGD  context  are  [7,8,10,11,12]. 
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§2.  Fundamentals  of  Laguerre  Geometry 

In  the  cyclographic  model  of  Laguerre  geometry,  a sphere  S with  center 
m — (mi, m2,  m3)  and  radius  r is  described  as  a point  s = (mi,  m2,  ms,?-)^  £ 
H ’ . It  is  an  intuitive  model,  since  the  space  of  points  (spheres  with  vanishing 
radius)  is  embedded  in  the  cyclographic  model  as  the  hyperplane  X4  = 0.  The 
absolute  quadric  fl,  which  reads  in  homogeneous  coordinates  {xjx^  :=  Xi)  as 
fl  : £0  = X 1 + X2  + £3  — £4  = 0 plays  an  important  role.  It  defines  a pseudo 
Euclidean  (pe)  metric  in  via  the  scalar  product 


(fl)b)pe  ■—  + 03^*3  — 04^4’  (3) 

The  pe  distance  ||a  — 6||pe  :=  y/{a  — h,a  — b)pe  of  two  points  o,  6 € R®’^ 
measures  the  tangential  distance  of  two  spheres  A,B.  First  we  study  the 
most  simple  family  of  spheres  described  as  a line  L with  direction  vector  a 
in  R®’^.  There  are  three  Ccises:  If  {a,a)pg  > 0,  then  L is  called  a hyperbolic 
line  corresponding  to  a family  of  spheres  whose  envelope  surface  is  a right 
circular  cone  (Fig.  2,  left).  If  {a,a)pe  = 0,  then  L is  called  a parabolic  or 
isotropic  line  belonging  to  spheres  in  oriented  contact.  If  (a,  o)pe  < 0,  then  L 
is  called  a elliptic  line  and  two  arbitrary  spheres  of  such  a family  do  not  have 
a common  tangent  plane.  These  three  cases  of  sphere  families  are  plotted 
in  Fig.  1 (but  to  simplify  matters  in  R^  instead  of  R^).  A canal  surface 
is  completely  determined  by  the  set  of  its  tangent  cones.  They  belong  to 
hyperbolic  tangent  vectors  of  a curve  c(t)  e R®’^  Vice  versa,  the  set  of 
hyperbolic  tangent  vectors  of  c{t)  corresponds  to  the  set  of  tangent  cones  of  a 
canal  surface  along  the  characteristic  circles.  For  elliptic  tangent  vectors,  there 
does  not  exist  a real  tangent  cone.  Therefore,  that  case  has  to  be  avoided  in 
the  specification  of  canal  surfaces.  Discrete  parameter  values  to  with  parabolic 
tangent  vectors  (||c(to)||pe  = 0)  are  tolerated.  At  these  parameter  values,  the 
canal  surface  is  closed  (Fig.  2 right). 

Corollary  1.  Any  real  -continuous  canal  surface  C in  R^  can  be  described 
as  a G^-continuous  curve  c(t)  £ R^’^  with  hyperbolic  tangent  vectors. 

A Laguerre  transformation  r in  R^'^  is  an  affine  mapping  x XAx  + b, 
0 < A e R which  preserves  fl.  It  turns  out  that  this  condition  is  fulfilled  iff 
A is  an  orthogonal  matrix  with  respect  to  the  pe  metric,  i.e. 

A'^EpeA  = Epe  :=  diag{\,  1, 1,  -1). 


(4) 
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Fig.  2.  Left:  Tangent  cone  of  two  spheres;  right:  Effect  of  parabolic  tangent  vector. 

Therefore,  Laguerre  transformations  are  pe  similarities  in  Since  the 

Euclidean  space  is  embedded  in  R®’^,  any  Euclidean  translation  or  rotation 
is  a Laguerre  transformation.  The  offset  operation  with  offset  distance  d is 
described  as  a translation  rj  in  the  o:4-direction  (A  = 1,  .4  = E,b  = (0, 0, 0,  d)). 
A further  example  is  the  pe  rotation  around  a fixed  2-plane.  If  the  X2,  xs-plane 
is  fixed,  the  pe  rotation  is  represented  by 

. f xi\  ( cosha  sinha'\  f xi\ 

J ysinha  cosha  J ^0:4  J ' ' ^ 

Any  Laguerre  transformation  with  A = 1 preserves  the  tangential  distance. 
The  following  examples  show  how  Laguerre  transformations  act  on  simple 
curves  of  R^’^. 

Example  1.  Consider  a straight  line  L € R^.  The  image  line  under  the  offset 
mapping  describes  a cylinder.  An  additional  pe  rotation  Tr  oTd(L)  yields  a 
hyperbolic  line  corresponding  to  a right  circular  cone.  Vice  versa,  any  hyper- 
bolic line  in  R^’^  can  be  mapped  via  an  appropriate  Laguerre  transformation 
to  L. 

The  line  in  Example  1 does  not  change  its  type  under  the  mappings  Tr  and 
Td.  Since  the  sign  of  (a,  o)pe  is  not  modified  by  any  Laguerre  transformation, 
we  obtain: 

Corollary  2.  The  type  (hyperbolic,  parabolic,  elliptic)  of  straight  lines  in 
the  cyclographic  model  R^’^  is  invariant  under  Laguerre  transformations. 

Example  2.  Consider  a Euclidean  circle  K S R^.  Its  image  Td{K)  describes 
a torus,  and  the  pe  circle  Tr  ° Td{K)  corresponds  to  a Dupin  cyclide.  Any 
Dupin  cyclide  can  be  interpreted  as  a circle  in  R^’^  with  respect  to  the  pseudo 
Euclidean  metric.  For  more  details  and  a proof  of  this  fact,  see  [7]. 

The  pe  circles  might  be  utilized  to  describe  G^-continuity  between  canal 
surfaces:  The  osculating  pe  circle  of  a curve  c(t)  S R^’^  belongs  to  the  oscu- 
lating cyclide  of  a canal  surface  G,  which  specifies  the  curvature  behaviour  of 
G.  A technique  to  compute  osculating  circles  of  space  curves  in  the  Euclidean 
and  non- Euclidean  space  is  given  in  [1]. 
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§3.  Isotropic  Hypersurfaces 

Definition  1.  Let  C C be  an  oriented  -continuous  canal  surface,  which 
belongs  to  the  curve  c e R^’^.  The  isotropic  hypersurface  r(c)  C R®’^  is  the 
union  of  all  points  in  R^’^  corresponding  to  all  oriented  tangent  spheres  of  C. 

For  a fixed  surface  point  of  C,  all  oriented  tangent  spheres  are  represented 
as  isotropic  lines  in  R^’'^.  Therefore,  F(c)  is  formed  by  a two-parametric 
family  of  isotropic  lines.  Applying  Corollary  2 to  the  definition  of  isotropic 
hypersurfaces,  we  obtain: 

Proposition  1.  Isotropic  hypersurfaces  are  invariant  under  Laguerre  trans- 
formations T(r(c))  = r(r(c)). 

The  simplest  example  is  the  isotropic  hypersurface  of  a sphere  5 specified 
by  a point  s € R^’^.  It  is  a hypercone  T(s)  = {x  — s,x  — s)pc:  = 0.  Now 
consider  a canal  surface  C characterized  by  a curve  c(t)  G R^’^  with  hyperbolic 
tangent  vector  c{t).  As  shown  in  [7],  its  isotropic  hypersurface  can  be  achieved 
in  a two-step  procedure:  First  intersect  the  pe  polar  plane  Cpg(t)  with  the 
absolute  quadric  0,  (the  pe  polar  plane  of  a vector  o € R^’^  is  defined  as 
Opg  : oTEpeX  = 0).  Then  join  the  resulting  family  of  conics  k{t,s)  and  the 
curve  points  c(<)  with  straight  lines.  Or  as  a formula: 

F(c)  =U(  c(t)*{c^^(t)nn}, 

where  A-k  B denotes  all  straight  lines  joining  A and  B. 

Isotropic  hypersurfaces  have  several  properties  which  can  be  applied  in 
the  context  of  geometric  design: 

Proposition  2. 

1)  The  intersection  of  the  isotropic  hypersurface  F(c)  with  R^  generates  the 
canal  surface  C itself:  F(c)  fl  R^  = C. 

2)  The  intersection  of  rd(F(c))  with  R^  gives  the  offset  surface  of  C. 

3)  Consider  two  canal  surfaces  C\,C2  € R^  with  corresponding  curves 
Ci(t),C2(t)  € R^’^.  These  canal  surfaces  touch  each  other  if  C2  is  lying  on 
F(ci)  (resp.  if  Ci  is  lying  on  F(c2)j. 

4)  The  bisector  surface  of  Cj  and  C2  is  formed  by  the  orthogonal  projection 
of{F(ci)nr(c2)}  toR^ 

Proof:  1)  The  tangent  spheres  with  radius  zero  of  a canal  surface  are  pre- 
cisely its  surface  points.  They  are  contained  in  the  hyperplane  R^  of  the 
cyclographic  model.  2)  follows  from  the  offset  property  of  Td  and  Proposition 
1.  3)  follows  immediately  from  Definition  1.  To  prove  4),  note  that  the  bisec- 
tor is  characterized  as  the  center  of  spheres  which  have  oriented  contact  to  Ci 
and  C2-  In  R^’'^  these  spheres  are  characterized  by  {F(ci)  fl  F(c2)}.  □ 

We  complete  this  section  with  few  examples  which  make  use  of  Laguerre 
transformations  and  isotropic  hypersurfaces. 
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Fig.  3.  Bisector  of  2 spheres,  resp.  of  sphere  and  canal  surface. 

Example  3.  Using  the  previous  proposition  and  Example  2,  implicit  and 
parametric  representation  of  cyclides  C can  be  derived  easily.  With  a Laguerre 
transformation  r,  the  corresponding  pe  circle  can  be  mapped  to  a simple 
normal  form  ^=^  = Xi+X2~ ^ = 0 of  a image  cyclide  C.  Its  isotropic 
bypersurface  reads  in  implicit  form  as 

r(c)  : (*o  + - 4xg(xi  + xl)  = 0, 

and  in  parametric  form  as 

(Xq(u,  V,  w),Xi(u,  V,  w),X2(u,  V,  w),  Xg(u,  V,  w),  X^{u,  V,  w))  = 

([1  - u;][l  + 1)^]  + w[l  - u^],  [1  - w][l  - v\  2[1  - w]v,  2wu,  - w[l  + u^]). 

The  implicit  and  parametric  representation  of  the  original  cyclide  C can  be 
obtained  from  {r“^(r(c))}nR^.  The  intersection  procedure  is  straightforward, 
since  r is  a linear  mapping.  In  the  parametric  case,  the  parameter  w can  be 
eliminated  easily,  because  it  occurs  linearly  in  the  parameter  representation. 

Example  4.  Consider  two  spheres  81,82  specified  by  points  si,«2  € 

The  intersection  r(si)  fl  F(52)  of  their  isotropic  hypersurfaces  is  contained  in 
the  hyperplane  H^2  ' {x  — (si  + S2)/2,  S2  — Si)pe  = 0-  Thus,  the  bisector  of  81 
and  82  is  the  orthogonal  projection  of  H12  D r(si)  onto  R^.  Since  r(si)  is  a 
hypercone,  the  bisector  surface  is  a quadric  surface.  Fig.  3 shows  two  spheres 
and  a part  of  their  bisector  surface  (hyperboloid  of  two  sheets). 

Example  5.  A similar  technique  allows  the  computation  of  the  bisector  of  a 
sphere  8 and  a canal  surface  C (speciGed  by  s and  c(t)  € Now,  Hi2(t)  : 

{x  — {s  + c(t))/2,  c{f)  — s)pe  = 0 depends  on  t.  The  bisector  surface  is  obtained 
by  the  orthogonal  projection  of  H12  D r(c)  onto  R^.  For  a rational  canal 
surface,  F(c)  is  a rational  ruled  hypersurface  with  parameter  representation 
g{t,u,v)  in  R^’^  which  is  linear  in  v.  After  the  intersection  and  projection 
operation,  we  obtain  a rational  bisector  surface  b{t,u).  Fig.  3 depicts  the 
bisector  b{t,  u)  of  a sphere  and  a canal  surface  c(t,  u)  with  cubic  spine  curve 
and  rational  parametric  degree  (7,2).  b(t,u)  has  parametric  degree  (10,2). 

Rational  bisectors  of  spheres  and  PN  surfaces  are  studied  in  [3].  Fur- 
ther applications  of  the  concept  of  isotropic  hypersurfaces  can  be  found  in  [8] 
(tangent  plane  property)  and  [12]  (offset  property). 
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Fig.  4.  Cyclide  and  cubic  canal  surface  with  control  spheres, 

§4.  Control  Sphere  Representation 
Consider  a Bezier  curve 

= (6) 

k=0 

of  degree  n in  the  cyclographic  model  with  control  points  sjt  6 They 

correspond  to  control  spheres  of  the  canal  surface  C described  by  x{t).  Well- 
known  properties  of  the  Bernstein-Bezier  technique  can  be  applied  directly  to 
the  Laguerre  geometric  approach  (this  idea  was  introduced  in  [12]): 

Proposition  3. 

1)  The  de  Casteljau  algorithm  enables  a stable  computation  of  the  moving 
spheres  of  a canal  surface  and  can  be  used  to  subdivide  it  into  two  parts. 

2)  The  control  spheres  So  and  Si  (S„  and  S„-i)  generate  the  tangent  cone 
of  the  canal  surface  at  t = 0 (t  = 1). 

3)  A (rational)  canal  surface  (with  positive  weights)  lies  in  the  convex  hull 
of  its  control  spheres. 

The  control  structure  is  useful  from  the  designers  point  of  view  because 
the  influence  of  moving  the  control  spheres  or  changing  their  radius  or  weights 
is  analogous  to  the  well-known  curve  case.  For  example,  Fig.  5 views  the 
influence  of  modifying  one  weight  of  a canal  surface  with  cubic  spine  curve: 
The  weight  /?i  is  increased  from  1 (left)  to  10  (right).  All  other  weights  satisfy 

A = l. 

We  have  seen  some  advantages  of  using  the  control  sphere  representation 
of  canal  surfaces.  However,  often  it  is  necessary  to  know  a (rational)  tensor- 
product  representation.  Therefore,  one  has  to  analyze  the  correlation  between 
control  spheres  and  classical  control  points  of  a canal  surface.  For  the  cyclide 
case  there  exists  a simple  geometric  relation,  which  is  described  in  [8].  In 
the  case  of  an  arbitrary  canal  surface,  the  problem  was  solved  by  Pottmann 
and  Peternell  [11],  They  proved  the  surprising  result  that  any  canal  surface 
with  rational  spine  curve  m{t)  and  rational  radius  function  r{t)  is  rational. 
Furthermore,  they  proved  that  the  problem  of  finding  the  rational  tensor- 
product  representation  can  be  reduced  to  the  problem  of  finding  two  rational 
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Fig.  5.  Influence  of  different  weights. 

functions  pi{t)  and  p2{t)  which  satisfy 

pI  + pI=  p:=ml+ml  + ml-r^.  (7) 

Although  the  existence  of  solutions  of  (7)  can  be  proved  via  factorization  over 
the  complex  field,  it  is  not  trivial  to  find  the  (non-unique)  solutions.  There 
are  two  different  options  to  perform  the  conversion  from  control  spheres  to 
control  points: 

1)  Use  arbitrary  rational  curves  x{t)  € to  describe  canal  surfaces.  Then 
it  is  impossible  to  solve  (7)  exactly  and  numerical  methods  are  required 
to  compute  the  control  points  fi:om  the  control  spheres  of  the  dedicated 
canal  surface. 

2)  Use  rational  curves  x{t)  € which  have  the  property  that  p{t)  is 
a square  in  the  polynomial  ring  R[t].  Then  an  exact  conversion  from 
control  spheres  to  control  points  can  be  realized  using  the  algorithm  of 
Pottmann  and  Peternell  [11]. 

In  this  paper  we  will  follow  the  second  option.  It  results  in  the  concept 
of  Minkowski  Pythagorean  Hodograph  curves. 


§5.  Minkowski  Pythagorean  Hodograph  Curves  in  R^’^ 

Planar  and  spatial  Pythagorean  hodograph  curves  have  been  introduced  by 
Farouki  and  Sakkalis  [4,5].  Recently  Moon  [9]  has  generalized  this  class  of 
curves  and  investigated  the  Minkowski  Pythagorean  hodograph  (MPH)  curves 
of  R^’^.  Here  we  need  a further  generalization:  the  MPH  curves  of  R^’^: 

Definition  2.  A polynomial  (rational)  curve  x{t)  = (xi(t),  X2{t),  X3(t),  X4(t)) 
in  R^’'^  is  a Minkowski  Pythagorean  hodograph  (MPH)  curve  if  there  exists  a 
real  polynomial  (rational  function)  a(t)  which  satisfies 

s?(t)  + x^(t)  + xl(t)  - x\(t)  = a^(t).  (8) 

As  described  before,  for  canal  surfaces  described  via  MPH  curves  of  R^’^, 
an  exact  conversion  from  control  spheres  to  control  points  can  be  realized.  We 
now  present  some  additional  important  properties  of  MPH  curves. 
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Proposition  4. 

1)  The  tangent  vector  of  a MPH  curve  is  hyperbolic  except  a finite  number  of 
parabolic  tangent  vectors,  i.e.,  the  corresponding  canal  surface  is  always 
real. 

2)  The  MPH  property  is  invariant  under  Laguerre  transformations. 

3)  PH  space  curves  are  exactly  those  MPH  curves,  which  hold  X4  = 0 

Pe  circles  are  examples  for  rational  MPH  curves  of  degree  2.  In  order  to 
create  polynomial  MPH  curves,  equation  (8)  could  be  solved  in  a polynomial 
ring  ]R[t].  A possible  solution  (determined  via  stereographic  projection)  is 

(il(t),  X2{t),  X3{t),  Xi{t), 

{2uqUi,  2uqU2,  ~ul  + uf  + Jij  “ W3,  2uoU3,  Mq  + wf  + — ti|). 

The  polynomials  u{t)  :=  [uo(t),ui{t),U2{t),U3{t))  of  degree  n are  mapped  via 
(9)  to  polynomials  of  degree  2n.  Taking  account  of  Xi{t)  = f Xi(t)  dt  + C 
with  i = 1, . . . , 4,  we  obtain  a polynomial  MPH  curve  of  degree  2n  + 1.  Due 
to  space  limitations,  we  postpone  a more  detailed  description.  However,  the 
basic  principles  of  the  analytic  construction  are  the  same  as  for  MPH  curves 
of  which  are  inspected  by  Moon  [9]  in  detail. 

An  alternative  approach  is  to  construct  MPH  curves  as  Laguerre  images 
of  PH  curves.  For  example,  one  could  make  use  of  well-known  properties 
on  spatial  PH  cubics  [5,6,14]  for  building  MPH  cubics.  Consider  a cubic 
MPH  Bezier  curve  x(t)  with  control  points  6o.5ii 52,5s  S 5q  and  5i 

resp.  52  and  5s  define  two  hyperbolic  tangent  vectors  ti  and  <2,  which  span  a 
three-dimensional  plane  H^.  If  there  exits  a Laguerre  transformation  r with 
r(H^)  = R^,  then  x(t)  = r~^(y(t))  is  a Laguerre  image  of  a (spatial)  PH 
cubic  y(t)  € R^  (that  fact  comes  from  properties  2)  and  3)  of  proposition  4 
in  combination  with  the  convex  hull  property  of  Bezier  curves). 

Proposition  5.  Consider  a hyperplane  H^  C R®’*  spanned  by  two  skew 
hyperbolic  lines  ti  and  t2-  There  exists  a Laguerre  transformation  t with 
r(H^)  = R^  iff  all  lines  joining  ti  and  <2  are  hyperbolic. 

Proof:  All  straight  lines  in  R^  are  hyperbolic.  Because  of  Corollary  2,  it 
is  necessary  for  the  existence  of  r,  that  any  line  of  the  linear  congruence 
ti  -k  t2  is  hyperbolic,  too.  The  sufficiency  is  proved  by  construction:  With  a 
pe  rotation  ti  one  hyperbolic  line  ti  can  be  mapped  to  R^.  Ti(ti)  and  the 
point  p :=  Ti{t2)  n R^  span  a 2-dimensional  plane  H^  C R®.  An  appropriate 
Euclidean  motion  T2  maps  H^  onto  the  X2,X3-plane.  Let  q = (91,0,0,94) 
be  the  intersection  point  of  T2  ori(t2)  with  the  Xi,X4-plane.  Case  1)  |9i|  < 
I94I:  the  pe  rotation  (5)  with  a =arctanh(— 91/94)  maps  T2  o ri(H^)  onto 
the  hyperplane  xi  = 0.  Since  it  contains  elliptic  lines,  this  contradicts  the 
assumption  that  ti  * <2  is  hyperbolic.  Case  2)  |9i|  = I94I  : T2  o ti(H®)  is 
the  hyperplane  xi  = X4.  Since  it  contains  parabolic  lines,  we  also  get  a 
contradiction.  It  remains  to  consider  Case  3)  [91 1 > I94I:  the  pe  rotation  (5) 
with  a =arctanh(— 94/91)  maps  T2  oti(H^)  onto  R^.  □ 
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Fig.  6.  Interpolation  of  spheres  with  cyclide  spline. 

§6.  Interpolation  with  Canal  Surfaces 

Finally  we  briefly  discuss  simple  interpolation  problems  concerning  canal  sur- 
faces. The  construction  of  canal  surfaces  which  interpolate  given  spheres  or 
tangent  cones  is  reduced  in  the  cyclographic  model  to  finding  an  interpolating 
curve  in  IR^’''.  Therefore,  standard  algorithms  can  be  generalized  from  E.® 
to  in  order  to  solve  the  problem.  Indeed,  one  has  to  take  into  account 
that  these  methods  could  generate  curves  with  elliptic  tangent  vectors,  which 
does  not  comply  with  real  canal  surfaces!  This  problem  can  be  avoided  using 
algorithms  which  deal  with  MPH  curves. 

Example  6.  n -I- 1 spheres  P,  (i  = 0, . . . ,n)  can  be  interpolated  with  a G^- 
continuous  cyclide  spline  (Fig.  6).  Each  piece  is  characterized  by  a pe  circle 
segment  with  control  points  bk,i  6 E^’^  and  weights  pk,i  € E,  fc  = 0, 1, 2. 
Per  given  tangent  vector  at  the  segment  is  determined  uniquely. 
Thus  with  any  starting  vector  foi  the  pe  circular  spline  curve  can  be  produced 
successively.  To  perform  the  computation,  one  can  generalize  an  ordinary 
circular  spline  algorithm.  The  formulas  (10)  are  cited  from  [13],  just  replacing 
the  Euclidean  metric  by  the  pe  metric: 


^0,1 


= Pi-1 1 Ki 

Po,i  = 1, 


, IIAPill^e 

_ (^Pii^i-l)pe 
ll^Pillpell^i-llIpe 


^2,i  — Pi, 
P2,i  = 1. 


(10) 


The  only  condition  on  the  input  data  is  l|Apj||pe  :=  ||pj  — Pi_i||  > 0,  i.e.,  the 
spheres  Pj_i  and  Pj  cannot  lie  inside  each  other. 

Example  7.  A piecewise  Hermite  interpolant  can  be  realized  with  a cu- 
bic MPH  canal  surface.  Consider  two  spheres  Pj_i,  P,  (specified  by  points 
Pi-1)  Pi  ill  E^’'^)  and  tangent  cones  Tj_i,  Tj  (described  as  hyperbolic  lines 
g E^’^  spanning  a hyperplane  C E^’^).  If  this  input  data  satisfy 
the  condition  of  Proposition  5,  then  there  exists  a Laguerre  transformation 
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T which  maps  11^  to  and  a well-known  Hermite  interpolation  algorithm 
[6,14]  produces  a cubic  PH  Hermite  interpolant  c(<)  £ IR^.  Its  Laguerre  image 
T”'^(c(t))  is  a cubic  MPH  curve  and  interpolates  the  original  input  data. 
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